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The NA61/SHINE experiment aims to discover the critical point of 
strongly interacting matter and study properties of the onset of deconfinc- 
ment. It also performs precise hadron production measurements for the 
neutrino and cosmic rays experiments. These goals are being achieved by 
measurements of hadron production properties in nucleus-nucleus, proton- 
proton and proton/pion- nucleus interactions as a function of collision en- 
ergy and size of the colliding nuclei. This contribution presents preliminary 
results from the NA61 ion program on single-particle spectra and identi- 
fied particle multiplicity fluctuations in p+p interactions at the CERN SPS. 
Comparisons with results from p+p and Pb+Pb collisions obtained by the 
NA49 experiment are shown. 

PACS numbers: 25.75.Ag, 25.75.Nq, 25.75Gz 

1. Introduction 

The NA61 /SHINE experiment and its predecessor NA49 at the CERN 
SPS study an important region of the phase diagram of strongly interact- 
ing matter. The Statistical Model of the Early Stage (SMES) of nucleus- 
nucleus collisions [1J predicted the energy threshold for deconfinement at 
low SPS energies. Several structures in the excitation functions were ex- 
pected within the SMES: a kink in the pion yield per participant nucleon 
(change of slope due to increased entropy production as a consequence of 
the activation of partonic degrees of freedom), a sharp peak (horn) in the 
strangeness to entropy ratio, and a step in the inverse slope parameter of 
transverse mass spectra (constant temperature and pressure in a mixed 
phase). These signatures were observed in central Pb+Pb collisions by 
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NA49 around y/sNN = 7.6 GeV [2]. Fluctuation analysis may provide ad- 
ditional evidence of the onset of deconfinement. 

Moreover, theoretical considerations suggest a critical point (CP) of strongly 
interacting matter which may be observable in the SPS energy range accord- 
ing to most lattice QCD calculations [3]. Fluctuations and correlations are 
basic tools to study this phenomenon. Enlarged fluctuations are expected 
close to the critical point. In nucleus- nucleus collisions a maximum of fluc- 
tuations is expected when freeze-out happens near the CP. 

2. The NA61/SHINE experiment 

NA61/SHINE is a fixed target experiment located in the North Area 
of CERN. SHINE stands for SPS Heavy Ion and Neutrino Experiment. 
It was approved in 2007. NA61 greatly profits from the long development 
of the CERN particle sources and the accelerator chain as well as the H2 
beam line of the CERN North Area. The latter has recently been modified 
to also serve as a fragment separator as needed to produce the Be beam for 
NA61. NA61/SHINE inherited the basic components from its predecessor, 
the NA49 experiment, and added important upgrades. The NA61 setup is 
presented on Fig. [I] (left). 




Fig. 1. The NA61 setup (left) and the energy - system size scan of NA61 in progress 
(right). Blue rectangles indicate data collected by NA49, circles indicate planned 
NA61 data, green circles are data points already recorded by NA61. 



The detector consists of five Time-Projection Chambers. Three of them 
(VTPC-1, VTPC-2 and GAP-TPC) are located in a magnetic field. Two 
others (MTPC-L and MTPC-R) are located downstream of the VTPCs. 
The TPCs are supplemented by three Time-of-Flight walls located behind 
the MTPCs. This setup allows for good particle identification and momen- 
tum resolution. 

The NA61/SHINE goal is to measure hadron production in p+p, p+A, h+A 
and A+A interactions at the SPS energies, in order to: 

• obtain precise data on hadron production (spectra) in p+C interac- 
tions at 31 GeV/c for the T2K experiment for precisely computing 
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the initial neutrino flux from the T2K target at J-PARC, 

• perform reference measurements of 7T+C and p+p interactions for 
cosmic-ray physics (Pierre-Auger and KASCADE experiments) for 
improving air shower simulations, 

• search for the critical point of strongly interacting matter, 

• study the properties of the onset of deconfinement, 

• study of high px particle production (energy dependence of Raa)- 

In order to study properties of the onset of deconfinement and search for the 
critical point NA61 performs for the first time in history a two dimensional 
scan in system size and energy. Already gathered systems and energies as 
well as planned measurements are presented in Fig. [T] (right) in the phase 
diagram of temperature and baryo-chemical potential. 

3. NA61/SHINE measurements for the onset of deconfinement 

and critical point 

The following results were obtained from p+p data collected in 2009 by 
NA61 /SHINE at = (6.^7.6 - 17.3 GeV. Event and track cuts were 

chosen to select only inelastic interactions and charged particles produced 
directly in the interactions or via strong and electromagnetic decays. 

3.1. NA61 results on single-particle spectra 

Results on single-particle spectra for ir~ mesons are obtained using two 
methods. First, by the h~ method which is based on the fact that the ma- 
jority of negatively charged particles are tt~ mesons. The contribution of 
other particles was estimated and subtracted using the VENUS and EPOS 
models. Second, by the dE/dx method which uses information on particle 
energy loss in the TPC gas in the relativistic rise region to identify particles. 
Results from both methods are corrected for feed-down from weak decays 
and detector effects using simulations. Out-of-target interactions are sub- 
tracted using events recorded with the empty liquid hydrogen target. 
Transverse mass spectra of tt~ in p+p interactions at 6.2 — 17.3 GeV are 
presented in Fig. [2] (left). For all energies rriT spectra are approximately 
exponential. The ratio of NA49 spectra for central Pb+Pb collisions [U [5] 
to those of NA61 for p+p interactions, both normalized to unity, is shown 
in Fig. [2] (right). The concave shape of the ratio seems to be energy inde- 
pendent. 

Rapidity spectra of 7r~ mesons in inelastic p+p interactions at 6.2 — 
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Fig. 2. Transverse mass spectra of ir mesons in inelastic p+p interactions at 
6.2 - 17.3 GeV (left) compared to central Pb+Pb data of NA49 (right). Statistical 
errors are indicated bv vertical bars and systematic uncertainties bv colored bands. 
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(it*) measured for NA49 
(7i-) = 60 + W for P 4.p 
(N w i mean number of wounded ni 
(Phys. Atom. Nucl. 64, 1841 (2001 )) 
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Fig. 3. Rapidity spectra of n~ mesons in inelastic p+p interactions at 6.2 — 17.3 
GeV (left) and mean pion multiplicity per wounded nucleon as a function of the 
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17.3 GeV are presented in Fig. [3] (left). An excellent fit at all energies 
is obtained by the sum of two symmetrically displaced Gaussian functions 
from which the mean tt~ multiplicity was derived. 

The mean yield of pions (tt + + 7r~ +7r°) per wounded nucleon (Fig. [3] (right)) 
was calculated using phenomenological isospin factors [6]. NA61/SHINE 
measurements are in agreement with the world data. The mean multiplic- 
ity in central Pb+Pb rises faster than in p+p collisions with a crossover 
at about ^snn = 10 GeV. The precision of the data allows for a detailed 
study of the onset of deconfinement. 

Rapidity spectra of tt~ , tt + , K~ mesons and protons produced in inelas- 
tic p+p interactions at 17.3 GeV are presented in Fig. |4j The spectra are 
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compared with the corresponding NA49 results [TJ El Ej^] as well as with 
the EPOS model [lOj . The experimental data agree for all particle types, 
between both analysis methods and NA61 and NA49. The EPOS model 
shows significant deviations from the data for K~ mesons and protons. 




y y 

Fig. 4. Rapidity spectra of 7r~, 7r + , K~ mesons and protons in inelastic p+p 
interactions at 17.3 GeV. The NA61 results are compared with the corresponding 
NA49 data and predictions of the EPOS model. 



3.2. NA61 results on chemical fluctuations 

Particle identification for the chemical fluctuation analysis is based on 
energy loss measurements in the relativistic rise region. It does not al- 
low unique identification of hadrons. The mean hadron multiplicities are 
obtained by fitting a sum of response functions for each particle type to 
the measured inclusive dE/dx distributions in narrow phase space bins. In 
order to analyze fluctuations of identified hadrons in p+p interactions of 
NA61/SHINE, as well as Pb+Pb collisions of NA49 a new approach, the 
so-called identity method [TTJ [T2J [13] , was used. It allows to obtain sec- 
ond and third moments (pure and mixed) of identified particle multiplicity 
distributions corrected for the effect of imprecise event-by-event particle 



2 NA49 points are obtained by interpolation from xf — Pt spectra 



6 



MMP'PROC PRINTED ON JANUARY 1, 2013 



identification. 

Two quantities, scaled variance Wj of the multiplicity distribution and the 




Fig. 5. Collision energy dependence of the scaled variance of identified hadrons 
multiplicity distributions in inelastic p+p interactions (left). The $^ fluctuation 
measure as a function of collision energy for p+p (NA61) and central Pb+Pb 
(NA49) collisions (right). The p+p data are compared with the EPOS model 
predictions. 



two-particle measure were chosen for the chemical fluctuation analysis. 
The scaled variance is defined as: 

u - < jg > ~ < N * > 2 m 

where < iVj > and < iVj > 2 are the mean multiplicity and second moment 
of the multiplicity distribution of particles of type i, respectively. 
The scaled variance is an intensive measure [H] i.e. it is independent of the 
number of wounded nucleons in the Wounded Nucleon Model or volume in 
the Grand Canonical Ensemble but it depends on their fluctuations. The 
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last feature makes it difficult to compare p+p interactions with nucleus- 
nucleus collisions. For the Poisson distribution uj = 1. 
Figure [5] (right) shows the scaled variance for it = ir + + vr~, K = K + + K~ 
and p = p + p as a function of collision energy. Results refer to the NA61 
acceptance for chemical fluctuation analysis |15| . 

The scaled variance for all studied hadron types increases with increasing 
collision energy. The trend as well as the magnitude of the effect is well 
reproduced by the EPOS model. As pions are the most abundant particles 
their fluctuations should be dominated by the KNO scaling |16j . It yields 
a linear increase of the scaled variance with mean multiplicity of charged 
particles. In NA61 this trend is modified by the limited detector acceptance. 
Values of ujk for all energies are close to or above one. This is most proba- 
bly caused by strangeness conservation which leads to a correlation between 
the production of K + and K~ mesons. Again the effect is weakened by the 
limited acceptance. To the contrary baryon number conservation seems 
to suppress proton fluctuations. This is because the proton multiplicity is 
mostly given by the two initial protons as production of proton-antiproton 
pairs is strongly suppressed by their large masses. This suppression effect 
is again weakened by the limited acceptance and it is expected to decrease 
with increasing energy. 

In order to compare results for p+p and central Pb+Pb collisions, the 
strongly intensive measure |144 [T7] defined for two hadron types, i and 
j, was chosen. It is defined as: 



where =< N { > ■ Uj + < Nj > ■ - 2 • (< Nj > - < iVj >< 
Nj >) • < Ni + Nj >. As a strongly intensive measure, <&ij is not only 
independent of number of wounded nucleons or volume but also of their 
fluctuations. 

Figure [5] shows the energy dependence of &ij for combinations of two hadron 
types: np, irK and pK. When no inter-particle correlations are present 
*y - 0. 

For ttK and pK increases with increasing energy. For $> np there is a 
minimum between 7.3 and 8.7 GeV. A similar but weaker effect is visible in 
Pb+Pb interactions. It also appears in the EPOS model. The increase of 
QttK for p+p interactions is not observed in Pb+Pb collisions. <3? p x shows a 
clear difference between results for p+p, which increase with increasing en- 
ergy, and Pb+Pb which decrease with increasing energy. Both dependencies 
cross zero at the same energy ^/snn ~ 8.7 GeV. 
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4. Conclusions 

The NA61 /SHINE experiment successfully started the planned 2D energy- 
system size scan with p+p interactions. The results are needed as reference 
for the onset of deconfinement study and the search for the critical point. 
The rriT spectra of tt~ mesons are approximately exponential. The com- 
parison to Pb+Pb interactions shows noticeable differences for the lowest 
and highest mx values which seem to be energy independent. A sum of two 
Gaussian functions is required to fit the rapidity spectra of 7r~ . 
The energy dependence of the mean pion multiplicity per wounded nucleon 
from NA61/SHINE agrees well with the previously established trend. 
Fluctuations of identified it, K and protons in inelastic p+p interactions are 
well described by the EPOS model. Conservation laws seem to play an im- 
portant role in fluctuations of hadron multiplicities. Comparison with data 
on central Pb+Pb collisions measured by NA49 shows large differences be- 
tween p+p and Pb+Pb in terms of fluctuations which are most pronounced 
for pK fluctuations. 
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